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Abstract 

We consider open XXX spins chain with two general boundary matrices submitted to one 
constraint, which is equivalent to the possibility to put the two matrices in a triangular form. We 
construct Bethe vectors from a generalized algebraic Bethe ansatz. As usual, the method also 
provides Bethe equations and transfer matrix eigenvalues. 

XXX spin chain [1] is one of the most studied integrable model. When the boundary conditions 
are the periodic ones, almost everything is known about it, starting from its spectrum up to 
asymptotics of correlation functions, and all the possible methods (e.g. coordinate [2] or algebraic 
Bethe ansatz [3, 4, 5, 6],...) have been successfully used. When the model has open boundary 
conditions, the situation changes drastically. Indeed, although the model is known to be integrable 
for general boundary matrices, for a long time only the case of diagonal boundary matrices was 
well- understood [7]. Recently [8, 9], a first step toward full resolution of the model was done, in the 
framework of algebraic or coordinate Bethe ansatz, and when one boundary matrix is triangular. 
Let us also mention [10] where the separation of variables has been used to tackle this problem. The 
aim of the present letter is to construct the algebraic Bethe ansatz for open XXX spin chain with 
two general boundary matrices with one relation among their entries. This relation is equivalent to 
the possibility to put the two matrices in a triangular form. 

When the two boundaries are triangular, the main difficulty lies in the fact that the number 
of excitations is not preserved anymore. To overcome this problem, we allow the ansatz to have 
states with unfixed number of excitations. The same idea was applied successfully in [9, 11, 12] in 
the context of coordinate Bethe ansatz. 

1 Reflection equation and transfer matrix 

The open XXX spin chain Hamiltonian can be constructed from the rational i?-matrix, satisfying 
the Yang-Baxter equation [13, 14], and if- matrices, satisfying the reflection equation [15, 7]. In this 
section, we will recall the main step of this construction and give the constraints on the parameters 
to get triangular boundary matrices. We also present its pseudo-vacuum state. 

'emails: samuel.belliard@univ-montp2.fr, Nicolas.Crampe@univ-montp2.fr, ragoucy@lapp.in2p3.fr. 
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1.1 Reflection equation 

The R-matrix associated to the XXX spin chain has the well-known form 

/ o(u) \ 

R(u) = b{u) ' 

c b(u) 
V a ( u ) J 

where 

a(u) = u + j] ; b(u) = u and c = rj . (2) 

This R-matrix, associated to the Yangian Y(sl 2 ), satisfies the Yang-Baxter equation [13, 14] 

R\ 2 (u - v)Ri 3 (u - w)R 23 (v - w) = R 23 (v - w)Ri 3 (u - w)Ri 2 (u - v) (3) 

where, as usual, the subscripts of the R-matrix indicate the spaces where the R-matrix acts non 
trivially. This R-matrix allows one to construct the monodromy matrix 

T a (u) = R al (u - a) Ra2(u - 6) ■ ■ ■ R aL {u - &) , (4) 

where L is the number of sites and £j are free parameters, called inhomogeneity parameters. This 
monodromy matrix is the cornerstone of the study of intcgrable periodic spin chains. 

To construct intcgrable spin chains with boundaries, we follow the method introduced in [15, 7] 
based on the reflection equations 



RM-u + v) B\{u) R 12 {-u -v-2ti) B* 2 (v) - B* 2 (v) Ri 2 (-u -v-2rj) b\{u) R 12 {-u + v) , (6) 



R 12 {u - v) B x {u) R 12 {u + v) B 2 {u) = B 2 {u) R 12 {u + v) B^u) R 12 {u - v) , (5) 

-% 

where . stands for the transposition. For later convenience, we introduce the following operators 

s^{u) = B 11 {u), #(u) = B ia (u) , tf(u) = fl ai (u) , @(u) = B 22 (u)--^—B 11 (u), (7) 

Zu + 1] 

where Bij(u) are the entries of the matrix B(u). The reflection equation (5) provides commutation 
relations between these operators. This computation is well-known [7] (see also [16]) and we report 
a list of these commutation relations in Appendix A. 

To construct an open spin chain, we need scalar solutions of these reflection equations. In the 
Y(sl 2 ) case, the most general scalar solutions are well-known and given, respectively for (5) and 
(6), by 

tsi \ { u(3 + a «7 \ . -jT, x / (—u — rj)f} + a (—u — rj)^ \ , . 

K ( u )=[ x a , and K ( u ) =[ / u I , ^ a i - ' ( 8 

y ' \ uo —up + a J V {—u — r])o (u + rjjp + a J 

where a,/3,7, S and a, /3, 7,(5 are free parameters. Using the monodromy matrix T(u) and the 
scalar solution K{u), we construct another solution of the reflection equation (5) via the dressing 
procedure 

B a (u) = T a (u) K a (u) T-\-u) . (9) 

We are now in position to introduce the transfer matrix associated to open XXX spin chain 

t{u) = tr a (K a (u) B a (u) ) , (10) 

which commutes for different spectral parameters (i.e. [t(u),t(v)] = 0). Finally, the intcgrable open 
XXX spin chain Hamiltonian is given by (in the case where £j = 0) 

H XXXopen = V ^lt(u)\ u = = g p,, +1 + ^(o) + -Lk< l (o). (ii) 
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1.2 Triangularization and constraints on the boundary parameters 

Although the spin chain with boundaries provided by the A"-matrices (8) is intcgrable (i.e. there 
exist L conserved charges), the computation of its eigenvectors remains an open problem. The 
main difficulty lies in the construction of a pseudo- vacuum, i.e. the determination of one particular 
Hamiltonian eigenvector. The case with diagonal boundaries does not share this difficulty and has 
been already treated in [7]. Then, using the i?-matrix invariance, the case where A and A can be 
simultaneously diagonalized has been treated in [17, 18]. The most general case treated up to now 
is when there exists a basis where A is triangular and A is diagonal [8]. In this paper, we propose 
a generalization of the algebraic Bethe ansatz (based on the ideas of papers [11, 9, 12]) to deal with 
the case when both A-matrices are triangular. 

Obviously [17, 18, 8], if we conjugate both A'-matrices by a constant matrix, the transfer matrix 
eigenvalues are unchanged. Then, we want to find a 2 by 2 matrix M such as M~ 1 K(u)M and 
M~ 1 K(u)M are upper triangular. Unfortunately, it is not always possible. It is a simple algebra 
exercise to show that one can do it if and only if the following constraint is valid 

(h-Sj) 2 -4((3j-^)(S/3~Sp)^0. (12) 

To our knowledge, it is the less restrictive constraint on the boundary parameters for which the 
eigenvalues and the eigenvectors of the transfer matrix are known (see Section 2.2). Evidently, we 
can deal with lower triangular matrices in the same way. 



1.3 Pseudo- vacuum 

From now on, we assume that the constraint (12) is satisfied. Then, we can triangularize the 
A-matrices (8) to get the following A-matrices 

ts( \ { ub + a uc \ , ~j7 1 \ ( -{u + rj)b + a -{u + rf)c. \ , . 

K M={ -ub + a) and K M={ (u + T,)b + a)> (13) 

where a, b, c and a, b, c are still free parameters. The relations between a, 6, c, a, b and c and the 
original parameters are given by 

a = a, b 2 =(3 2 + 1 5 1 c = 7 - [fi + b) (14) 

and similar relations with "bar" parameters. 

The transfer matrix to be diagonalized now reads 

t(u) = Ki(u) s?(u) + k 2 (u) 9{u) + n 12 (u) (15) 

where we used the notations 
2(u + rf) 

«i(w) = — (a — bu) , k 2 (u) = (u + rf) b + a , k%2 (u) = — (u + n) c. (16) 

2u + rj 

An important point is that for two triangular matrices, there still exists a simple eigenvector, 
called pseudo-vacuum. Indeed, let us consider the vector |f2) with L spin up i.e. 

i«>=(;)®(;)®-®(;) ■ ™ 

As explained in [8], when we choose A-matrices (13), this vector has the following properties 

£/(u)\n) = Ai(u) [O) with Aifa) = (a + bu)T\ ; (18) 

3 ii a (- u - 0) 

S(„)|n) = A 2 («) |0) with Aa(u) = Ma-Kv + v)) ^ Hu + QHu-^) 



2u + 1] a(u + £j)a(—u - 



)|fi) = 0. (20) 
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Using these properties, it is easy to show that |f2) is an eigenvector of t(u) 



t(u)\VL) = (ki(u)Ai(u) + k 2 (u)A 2 (u)) 



(21) 



2 Algebraic Bethe Ansatz 
2.1 Bethe vectors 

We are now in position to propose an ansatz for all the eigenvectors. Let us remark that if c 
vanishes, the ansatz used to study the case with both boundaries diagonal is still working. This 
trick was used previously in [8]. For the case with c different from 0, we need to generalize the 
ansatz. To this aim, we borrow the idea from the papers [11, 9, 12] where, instead of having a fixed 
number of excitations, we only fix N, the maximum of number of excitations. 

Before giving the explicit form of the Bethe eigenvectors, we need some definitions: 



rN 



X N = 



{ {h, t 2 , . . . , U} | 1 < ii < ■ ■ ■ < n < N } , 0<£<N 



N 



Us 



N 
I ) 



1=0 



{1,2,...,7V}\7, 



N 



(22) 
(23) 
(24) 

N 



Starting from a set of TV elements u = {u\, . . . , u^}, wc define uj = {ui \ i G /} for any I S X 
We can now define the Bethe vectors which depend on the Bethe parameters u = {u\, u 2 , ■ ■ ■ ujy}: 



i$ w (u)) = ^( u ) i B (^)>- 



lei" 



Wc have used the following vectors: 

|B(u/)> =Y[a(yn)\Sl) 



lei 



N 



(25) 



(26) 



iei 



and the following functions 



W/(u) 



n 



N 



A 2 (ui) 



c(2uj + rj) 
2{bui - a) 



J h(ui,u k ) 



I 6 1 



N 



(27) 



where the functions / and h arc defined in (59). 

Before giving the main result of this paper, let us comment the form of the Bethe vectors to 
clarify notations and conventions. We use the usual convention that a product over an empty set 
is equal to one, e.g. if I £ X$, we get 7 = and Wi = 1. We deduce that, for c — 0, we get 



JV 



(28) 



that corresponds to the usual ansatz, as expected, since for c — the left boundary becomes 
diagonal and the operator ^ is not present anymore in the transfer matrix (15). For N = 0, the 
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Bcthc vector |$°) reduces to the pseudo- vacuum |f2). Finally, let us give more explicitly the Bethc 
vectors for the cases N = 1,2: 

|$V)> = (#(ui) + A a (uO l { n 1 + 7 l\ ) |n> (29) 

\ 2(oui — a) J 

|* 2 (ui,u 2 )> = («(«i)«(«2) + W {2} (u u u 2 )^(u 2 ) + %(ui,ti 2 )%) + W % {u u u 2 )) | 0} (30) 



with 

Ti7 / s a / ,c(2ui+?7) («i -tt 2 +»7)(ui +U2 +2t?) , . 

W{ 2} (ui,u 2 ) = A 2 (wi) — = — — - — — — , (31) 

2{bu\ — a) {ui — u 2 ){ui +u 2 +rj) 

rxr / N A ( \ A I \ c( 2u l + V) c(2li 2 + 7/) [u X + U 2 + 2 7]) 

W0(ui,u 2 ) = A 2 (ui) A 2 (w 2 ) — = — — = — — - — - — r — . (32) 

2(oui — a) 2(bu 2 — a) [ui+u 2 ) 

2.2 Eigenvectors and eigenvalues 

We arrive at the main results of this paper. The Bethe vector (25) is an eigenvector of the transfer 
matrix t(u) defined by (15) i.e. we get 

t(u)\$ N (u)} = A(u)\$ N (u)) (33) 

with 

N N 

A(u) = ki(u)Ai(u) Y[f(u,u k ) + k 2 (u)A 2 {u) Y[h(u,u k ) (34) 

k=l k=l 

if the Bethe parameters {ui, u 2 , . . . , ujv} satisfy the following Bethe equations 

Al{Uk) -E(u k ){[lp^l (35) 



A 2 K) ^ k f{u k ,Uj) 



with 



2 u(a — bu) 

The proof of this result is given in Section 3. 

3 Proof of relation (33) 

3.1 Actions of si (u), @{u) and Sf(«) 

We need the actions of srf(u),@(u) and c €{u) on vectors of type |B(x)) = J[ £ i=1 38{xi)\£l) where 
x = {xi, . . . , Xi} can be any subset of Bethe parameters. The computation for the first two actions 
is a usual computation [7] using the commutation relations of Appendix A. We get 



*f(u)|B(x)) = A{(u,x)\E(x)) +J2 M k (u,x) 38(u) ]\SS( Xi )\n) (37) 

fc=i i=1 

l i 

#(u)|B(x)) = A|(«,x)|B(x)) + ^ N k (u,x) 3§{u) j[a{xi)\Sl) (38) 



fe=i <=i 
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where 

I i 

Af(«,x) = Ai(«) l[f(u,x k ) , A^,x)=A 2 (w) JJ/»(«,aj fc ) (39) 
fe=i fe=i 

and, using the notation x^j. = {a;!, . . . , Xfc_i, x^+i, ■ • ■ , x^} 

M k (u,x) = g(u,x k ) A e ~ 1 (x k ,x^ fc ) + w(u,x k ) A e 2 ~ 1 (x k ,x^ fc ) (40) 
N k (u,x) = k(u, x k ) A t f 1 (x k ,x^ fe ) + n(u,x k ) A[~ 1 (x k ,x^ fc ). (41) 

The computation of the action of ^(it) on |B(x)) is more involved, but after some algebra, using 
commutation relations of Appendix A, we get 

I i 
tf(u)|B(x)) = ]TG ? ( U ,x) |B(x^)) +^F ij (u,x) |B(x #iJ )>. (42) 

i— 1 2<j 

The functions F and G are given in terms of the functions defined in (60) and by the following 
relations: 

Gi(w,x) = A^ 1 (ii,x^,)^(m(w, Xi) + l(u,x t ))A^ 1 {x l ,y.^ i ) + p(u, Xj)A 2 _1 (x^x^jM 

+ A2 _1 (u,x #l )^((7(u,x t ) +y(u,x i ))A{~ 1 (x i ,x^ i ) + z(u, x l )A 2 ~ 1 (x l , x^H (43) 

i*y(u,x) = Ai" 2 (a; l ,x^ i j) (Zn(u, x i; x^A-f^fo-, x^) + Zi 2 (it, x { , Xj)A e 2 ~ 2 (xj, x^-)) 

+ A2 _2 (a; 4 ,x 7 i i j) (Zi 2 (w,x 3 ,2; l )Ai" 2 (x : ,-,x^ l . J ) + Z 22 (u, x i; 5G 3 -)A| _2 (a;j, x^)) (44) 

where 

7 , v Si] 2 x l Xj(x i + a;j)(M 2 - XjXj + 7?m) , . 

(ZXi + ??)(2Xj + J?) (X ? ; + + 77) (U + + 7]) (it + Xj + f])(U — Xi)(U — Xj) 

Z (u x- x ) - Wxj{x 3 - x,j + r])(u 2 + 7?m + XjXj + rjxj) 

(2xi + 77) (x, — Xj)(u + Xi + rf)(u + xj + T])(u — x,)(u — Xj) ' 

Z 22 (u x- x ) = 2ry 2 (^» + xj + 27/) (u 2 - (a, + r]){xj + 7/) + rju) 

(xi + Xj + rj)(u + Xi + rj){u + x 3 + rj)(u - Xi)(u - Xj) 

We want to prove relation (33) using relations (37), (38) and (42). We project the L.H.S and 
the R.H.S of (33) on all vectors |B(u/)), I E 1 N , and all vectors ^(u)|B(u/)), I £ 1 N , and prove 
that each relation we get in this way is true. 



3.2 Projection on |B(u)) 

This computation is similar to the usual one (with diagonal boundaries). Thus, we just sketch the 
proof. 

The only contributions of the L.H.S. of (33) are the first terms of relations (37) and (38) for 
£ = N. Using the explicit expressions (39), we recognize the R.H.S. of (33). This proves that the 
projection of (33) on |B(u)) is satisfied. 

3.3 Projection on &(u)\B(u^ k )) 

This computation is also similar to the usual one with diagonal boundaries. The only contributions 
of the L.H.S. of (33) are the second terms of the relations (37) and (38) (for I = N). There is no 
contribution from the R.H.S. of (33). Then, the relation to prove is 

Ki {u)M k (u,u) + n 2 {u)N k (u, u) = 0. (48) 

Using the explicit forms (40) and (41), the previous equation (48) is equivalent to the Bcthc 
equations (35). This proves the projection on <S?(u)|B(uj)) for / 6 3-n-v 
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3.4 Projection on |B(u 7 )) for / e If and I = 0, 1, . . . , N - 1 

This type of projection is not present in the usual diagonal boundary case. Therefore, we give here 
more details. There are contributions of the L.H.S. of (33) coming from the first terms of relations 
(37) and (38), from the first term of (42), and contributions of the R.H.S. of (33). The relation to 
prove is then 

W / (u)(k 1 (u)A £ 1 (m,u / ) + k 2 (w)A^(u,u / ) - A(u)) + K 12 {u)^2Wiuj(u) Gj(u, Ujuf) = 0. (49) 

In fact we are going to prove a stronger statement: the coefficient of A{(u, uj) in (49), as well as the 
coefficient of A^it, uj) in the same equation, both vanish identically. Indeed, the first coefficient 
reads 

coefi(u) = Ki(it)W/(u)(l - Y[f{u,Uj)) 

+Kl2(u) ^ W^uj(u)((to(u,-Uj) + l(u,Uj))A[(llj,Ul) + p(u,Uj)A2(%,u/)J (50) 

while the second one is just coef2(u) = cocfi(— u — r/). Thus, if one coefficient vanishes for all u, so 
does the other one. 

To show that cocfi(w) vanishes, we follow the technics used in [11, 9]: we prove that coefi(u) 
corresponds to the sum over all residues of some function. The function to consider is 

The poles of F (considered as a function of z) are located at z = uj, —uj — rj, 0, u and oo. It is a 
simple exercise to show that the sum over all the residues is then just equivalent to coefi(u) = 
when one replaces the functions Wj, A k , Ki, K12, m, I and p by their explicit expression (27), (39), 
(16) and (60), and uses the Bethe equations (35). This proves the projection on |B(u/)) for / £ Xf 
and 1 = 0,1,..., N - 1. 

3.5 Projection on ^(u)|B(u/)) for / G Xf and £ = 0, 1, . . . , N - 2 

This type of projection is also new in comparison to the case with diagonal boundaries. Therefore, 
we give also some details. There are contributions of the L.H.S. of (33) coming from the second 
terms of relations (37) and (38), from the second term of (42) and no contribution from the R.H.S. 
of (33). The relation to prove is then 

^V^ /U j(u)^ 1 (u)M i (u,u /ui ) + K 2 (")A r j(M,u /ui )^ +k 12 (u) ^ W IU{j , k} (u)F. lj (u,u IU{ . jM ) = 

jel j.fceT 
J j<k 

(52) 

The previous relation becomes the following functional relation when we used the explicit expres- 
sions of Wi, Mj, Nj, Fij Ki, K12 (see relations (27), (40), (41), (44) and (16)) and the Bethe 
equations (35) 

o = l(u, %) ( n ^ ~ n fa) + 2 s (^^' ~ uj ~ ??! ~ ui ~ ^ n f^f* + ( 53 ^ 

j(zj ke~ ke~ ke~ 

+Q(u,Uj,ut) Y\ h jk h lk + Q(u,-Uj - r),ui) fjkhk + Q(u,Uj,-Ui-r]) h jk fek) 

fceT keT k£J 

k^j,£ k^j,l k^j,£ 
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where hjk and fjk stand for h{uj,Uk) and f(uj,Uk) and 



L(u,Uj) = 2 (K 1 (u)g(u,u j ) + K 2 (u)n(u,u j )) E{ Uj ) — — ^ - (54) 

yzuj -t- r/ji^itj -t- ryj 



(6u J --o)(6« < -o)(ti J -+^ + 2j7) 
Q{v,,u h v,i) = -4 Zniu^^ui) ^(ujMue) —f- — — — — — -f — - — — . (55) 



To prove relation (53), we consider its R.H.S as a function of u, called X(u) and proves that it 
vanishes. Firstly, it is easy to see that X(u) is a rational function that tends to when u — > oo. 
Secondly, we remark that X(u) can possess poles only at u = Uj and u = —Uj — rj. Thirdly, we 
prove that its residues at u = Uj is equivalent to the sum over all the residues of the following 
function over z 

{z + Ujjiz-Uj 

and, by consequence, vanishes. We perform the same type of computation for the poles at it = 
— Uj — rj. This makes X{u) a rational function that vanishes at infinity and has no pole: it is equals 
to 0. This proves the projection on J?(w)|B(u/)) for I e and I = 0, 1, . . . , N — 2 and concludes 
the proof of relation (33). 



4 Conclusion 

We have constructed the algebraic Bethe ansatz for the XXX model with open boundary conditions 
characterized by two general boundary matrices with one constraint. This relation amounts to state 
that the two boundary matrices can be triangularized in the same basis. This property corresponds 
to the local gauge transformations used in XXZ model to diagonalize the boundaries, and possible 
only when constraints are applied to the boundary matrices. 

Several direction of investigations can follow: keeping the same model, one should compute the 
scalar products of (off-shell) Bethe vectors in order to have access to the correlation functions of 
the model; the same construction for the XXZ model with non-diagonal boundaries should also be 
done and compared with the previous works [19, 20, 21, 22]; generalization to different models such 
as spin chains based on algebras of higher rank or the Hubbard model should be also investigated; 
finally, the seek of an algebraic Bethe ansatz for general boundary matrices is a very exciting (but 
also very open) problem. 



A Commutation relations 

Using the reflection equation (5), we can find the exchange relations between the operators 2%, £ 
For our calculations, we only need the following ones 

\2${u),3g(v)\ = (57) 
af{u)S§{v) = f(u, v)S§{v)af{u) + g(u, v)^{u)^{v) + w(u, v)@(u)@(v), 

\€(u),£8(v)\ = m(u,v)£/(v)£/(u) + l(u,v)£/(u)£/(v) + q(u,v)j2/(v)S>(u) (58) 
+p(u, v)s?{u)2){v) + y(u, v)@(u)g/(v) + z(u, v)9(u)9(v) 
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with 



,/ \ (u - v - r))(u + v) 2nv -r\ 

f(u,v) = r, g(u,v) = -, w(u,v) = (59) 

(u + v + n)(u — v) (Zv + r])(u — v) (u + v + rj) 

(u - v + ?])(u + v + 2ij) -2rj(u + rj) 

h(u,v) = — — — , k(u,v) 



n(u, v) 



(u — v)(u + v + ?/) ' (u — v)(2u + rj) '' 

4 w rj (u + r)) 



(u + v + r))(2v + r])(2u + 7]) 



and 



f \ 2rju(u-v + ri) 2t?u 

m(u,v) = — — -, l(u,v) = — — — — (oO) 

v ' (2u + r))(u + v + r])(u-v) v ' (2u + rj)(2v + ri)(u-v) y ' 

T)(u + v) 2r)u 

1K U , V ) = 7 — ; ; — ^7 \> P( u > v ) 



(u + v + tj)(u — v) ' ' (2u + r])(u — v) 

y{u,V) = — ; ; — ; r, Z{U,V) ~ 



(u + v + rj) [2v + rj) u + v + rj 
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